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We develop a theory for interlayer pairing of chiral electrons in graphene materials which results 
in an unconventional superconducting (S) state with s-wave spin-triplet order parameter. In a 
pure bilayer graphene, this superconductivity exhibits a gap less property with an exotic effect of 
temperature-induced condensation causing an increase of the pairing amplitude (PA) with increasing 
temperature. We find that a finite doping opens a gap in the excitation spectrum and weakens 
this anomalous temperature-dependence. We further explore the possibility of realizing variety 
of pairing patterns with different topologies of the Fermi surface, by tuning the difference in the 
doping of the two layers. In trillayer graphene, the interlayer superconductivity is characterized 
by a two components order parameter which can be used to define two distinct phases in which 
only one of the components is non vanishing. For ABA stacking the stable state is determined by 
a competition between these two phases. By varying the relative amplitude of the corresponding 
coupling strenghes, a first order phase transition can occur between these two phases. For ABC 
stacking, we find that the two phases coexist with a possibility of a similar phase transition which 
turns out to be second order. 

PACS numbers: 74.78.-w, 74.70.Wz, 73.22.Pr 



I. INTRODUCTION 

The microscopic theory of superconductivity devel- 
oped by Bardeen, Cooper and Schrieffer 1 (BCS) is based 
on the pairing between electrons from opposite spin- 
subbands with the same Fermi surfaces. General- 
ization of the BCS pairing to composite systems of 
two (or more) types of fermions with different Fermi 
surfaces 2 has attracted great interest due to its appear- 
ance in various areas of physics including pairing of ul- 
tracold Fermi atoms^, color superconductivity in dense 
quark matter—, and neutron-proton pairing in nuclear 
systems 6 -. Despite the robustness of conventional metal- 
lic superconductivity^ against small mismatches between 
the paired Fermi surfaces, several unconventional pairing 
phases have been predicted to realize in the presence of 
larger mismatches. The most well known cases include 
the phase of breached pairing (BP) referring to a state 
of separated normal (N) and S phases in the momen- 
tum spaced, pairing between deformed Fermi surfaces 
with zero total momentum of Cooper pairs 9 and inho- 
mogeneous S phase of Larkin-Ovchinnikov-Fulde-Ferrel 
(LOFF)i£. These unconventional phases are resulted 
from different topologies of the Fermi surface. 

The search for new exotic S state is the frontier subject 
of low temperature physics, because it sheds light on the 
pairing symmetry and may lead to understanding pairing 
mechanism of superconductivity in high-T c n supercon- 
ductors and other newly discovered S materials^ 2 -. The 
recent discovery of graphene 13 i 14 , the two dimensional 
solid of carbon atoms with honeycomb lattice structure, 
and its associated bilayer and trilayer structures are ex- 
pected to provide still a new opportunity for realizing 
unconventional pairing states. Graphene has a specific 
zero-gap electronic band structure in which the charge 
carriers behave like 2D massless Dirac fermions with a 



pseudo-relativistic chiral property. In addition to reg- 
ular spin, electrons in graphene possess two additional 
quantum degrees of freedom, the so called pseudospin 
and valley. These features together with the fact that 
in graphene the carrier type, [electron-like (n) or hole- 
like (p)] and its density can be tuned conveniently, make 
these carbon-based material exceptional for realizing un- 
conventional superconductivity. The unusual features of 
superconductivity has been already predicted in mono- 
layer graphene 15 , where the pairing of electrons with op- 
posite sublattice pseudo-spin leads to appearance of an 
unusual spin-singlet p+ip-wave S phase with no gap in its 
excitation spectrum^. An intrinsic superconductivity, 
with plasmon or phonon mediated pairing interactions, 
can be realized in graphene coated with a metal 16,17 . For 
a coated bilayer and trilayer graphene, the formation of 
the S state is expected to be closely similar to that of 
the graphite intercalated with alkaline metals, for which 
a critical temperature up to 11.5 K has been reported^. 
Furthermore, progress has already been made in prox- 
imity induced superconductivity by fabrication of trans- 
parent contacts between a graphene monolayer and a su- 
perconductor (see for instance, RefsJ^). In a very re- 
cent study 20 , we have explored the exotic nature of in- 
terlayer superconductivity in pure bilayer graphene. It 
has been explained that how the interplay between the 
interlayer pairing of electrons with the same sublattice 
chirality and the asymmetric arrangement of the sublat- 
tices of the two layers result in a gapless superconductiv- 
ity with an unusual s-wave spin-triplet symmetry of the 
order parameter and anomalous thermodynamic proper- 
ties. We have obtained that the interlayer pairing allows 
for the possibility of a temperature-induced condensation 
causing an increase of the PA with increasing tempera- 
ture, and an entropy of the stable S state which can be 
higher than its value in N state. Motivated by these 
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findings for interlayer superconductivity of pure bilayer 
graphene, we study effect of a finite doping. Furthermore, 
considering the recent interest in the properties of tri- 
layer graphen o 21 i 22 , we extend our study to the interlayer 
pairing in trilayer graphene. For bilayer graphene, we 
find that the doping opens a gap in the excitation spec- 
trum which, in turn, weakens the temperature-induced 
condensation such that at high levels of the doping the 
temperature dependence of PA became similar to that 
of the conventional BCS gapi. We present phase dia- 
grams of interlayer superconductivity for both symmetri- 
cal and asymmetrical dopings with the two layers having 
the same and different levels of the doping, respectively. 

For trilayer graphene, we examine the interlayer pair- 
ing in two different types of ABA and ABC stacking in 
pure case and find an strong dependence on the type of 
the stacking. We show that the asymmetric ABC stack- 
ing can support stronger pairing gap than the symmetric 
ABA stacking. For ABA stacking, there is a competi- 
tion between two phases of interlayer superconductivity 
in which the pairing is realized between chiral electrons 
of the middle layer with those in only one of the other 
two layers. This depends on the relative values of the 
two corresponding coupling strengths. We explore the 
possibility of a phase transition between the two phases 
by varying the coupling strengths which turns out to be 
first order for ABA stacking but second order for ABC 
stacking. 

The paper is organized as follows. In the next section, 
we discuss the theoretical modeling of the interlayer su- 
perconductivity in graphene materials and study the re- 
sulted band structure with different topologies of Fermi 
surfaces. In Sec. we discuss the numerical results for 
the phase diagram and order parameter as function of the 
various parameters of both symmetrically and asymmet- 
rically doped bilayer graphene. Sec. [iVl is devoted to the 
case of interlayer superconductivity for trilayer graphene. 
Finally, in Sec. El we present our summary and conclu- 
sions. 



II. THEORY 

In order to study interlayer superconductivity, we con- 
sider a model based on graphene materials. We introduce 
the theory of interlayer superconductivity to the case of 
bilayer graphene and then this theory shall be developed 
to the case of trilayer graphene. Bilayer graphene is com- 
posed of two coupled graphene monolayers with A\ and 
B\ triangular sublattices in the top layer and A2 and B2 
triangular sublattices in the bottom layer according to 
Bernal stacking in which every A\ site of the top layer 
lies directly above a B2 site of the bottom layer. In 
the absence of superconductivity, the following Hamil- 
tonian can be used to describe tt electrons of the bilayer 
graphene^, 




FIG. 1: (color online) Left panel: Lattice structure of bi- 
layer graphene. Interlayer S correlations in the real space are 
shown by wavy lines which couple electrons of different lay- 
ers with opposite pseudospin degree of freedom. Right panel: 
Interlayer superconductivity in momentum space is shown by 
coupling of two time-reversed momentum states | — k) and \k) 
which are located in valleys K\ and K' 2 and belong to the dif- 
ferent layers. Since, interlayer pairing is antisymmetric with 
respect to the pseudospin degree of freedom, pairing between 
components \A\) and \B 2 ) is possible. 
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with the monolayer Hamiltonian, 

jjrnono = _ ^ ^ _ f („[, > ( . T + H .C.), (2) 

and the nearest neighbor interlayer hopping Hamiltonian, 
Hi B = -t ± J2(al ita b 2A , a + H.c), (3) 

where a\^^ (bi y i yCT ) and a\ ia (b\ i(T ) are the annihilation 
and the creation operators of an electron in the ith unit 
cell in the sublattice A(B) and the layer I a = zb denotes 
the spin state of the electron; n^ )Cr is the corresponding 
on-site particle density operator. The intralayer nearest 
neighbor hopping energy t ~ 3eV determines the Fermi 
velocity in graphene as ~ 10 6 m/s, and t± « OAeV 
(t±/t « 0.13) is the dominant interlayer hopping energy 
between the nearest neighbors A\ and B^\ the chemical 
potential \i\ can be controlled by a top- and bottom-gate 
electrodes independently. 

By introducing attractive interaction between the elec- 
trons of the sublattice A\ and B2 through the following 
potential, 

^ B = ~9±1 Yl a U* a hiM,i,a> b 2^>> ( 4 ) 
a, a' ,i 
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the interlayer superconductivity can be produced. Here, 
g±i is the S coupling energ y 17 ! 24 . The interaction poten- 
tial (|4]) has on-site local property in the 2D plane of the 
bilayer graphene. Therefore, the wave function of two- 
body problem associated with the two coupled electrons 
and the corresponding pairing potential are not extended 
in the bilayer plane. This results in the introduction 
of isotropic s-wave symmetry for the orbital part of the 
order parameter. Because the coupling occurs between 
the two electrons with the opposite pseudospin degree 
of freedom, Pauli exclusion principle imposes a pairing 
in the spin-triplet channel. Therefore, the total wave 
function, including the product of orbital, spin and pseu- 
dospin parts will be antisymmetric under the exchange 
of electrons. The following s-wave spin-triplet order pa- 
rameter can be used to decouple the total Hamiltonian 
of the bilayer graphene (H AB = H AB + Vf B ), 



(5) 



In the lattice space, as is shown in the left panel of Fig.HJ 
\,±i describes the interlayer pairing of the two electrons 
with the sublattices A1-B2, and the spins (and 
In the momentum space, for the case of t± =0 that 
the electronic band structures of the two monolayers are 
not affected by each other, there is a simple description 
for binding of the two time-reversed electronic states due 
to presence of attractive interaction as we discuss in the 
following. 

At low energies, the chiral momentum states \k) lo- 
cated in Dirac cone-shape band structures at the two 
inequivalent valleys are coherent superposition of sublat- 
tice pseudospin states \A) and \B), 



\k) = \A)+exp(i(p k )\B), 



(6) 



where cfk is the angel of the momentum direction. Inter- 
layer superconductivity induces partial pairing between 
the time-reversed momentum states \k) and \ — k) located 
in the Dirac cones of different layers such that the cou- 
pling only takes place between \Ai) and ^2) pseudospin 
parts of the electrons wave function in layers 1 and 2, 
respectively. Since the pairing of the time reversed chiral 
momentum states is partial, the interlayer superconduc- 
tivity (j5j) gives raise even parity of the order parameter 
in the frequency spaced. 

After decoupling of the interacting part of H AB by 
the order parameter ([5]) , one can diagonalize Hamiltonian 
in 2D momentum space. The general expression of the 
spectrum is large, but for t± = 0, we obtain a simple 
expression for the spectrum E^ which is, 
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FIG. 2: (color online) The plot of quasi-particle excitation 
spectrum E^ versus ek in S state. The dispersions correspond 
to the regimes of (a) BCS pairing with no effective Fermi 
surface for both branches / = ±1, (b) Sarma pairing with two 
effective Fermi surfaces for the branch I = —1 but no Fermi 
surface for the other branch I = 1, (c) pairing with single 
Fermi surface in I = —1 but no Fermi surface for the other 
branch 1=1, and (d) pn pairing with single Fermi surface in 
each branch / = ±1. 



h = (/ii — /42)/2 are free band dispersion, mean dop- 
ing and doping difference, respectively. Without loss of 
generality, we assume fii > fi2- Branch E^ t for any value 
of parameters h, /i, A_q, and I does not cross Fermi sur- 
face such that for breaking a paired state and sending 
quasiparticles to this state a finite energy is needed. This 
situation resembles metallic superconductor. Neverthe- 
less, E^ can cross Fermi energy through one point for 
I = 1 and either one or two points for I = — 1. As a re- 
sult, it generates gapless excitation modes. In Fig. [2] the 
quasiparticle excitation spectrum is shown. There are 
four possible features for E^ v Each feature with its own 
topology of Fermi surface can be accessible by tuning the 
value of h. The boundaries between four ranges of h asso- 
ciated with different topologies of the Fermi surfaces can 
be determined by analyzing the roots of E^ = which 
are three critical values, 




(9) 



A 2 ±1 + 2( M 2 + 4) + LA 



•44(Ai 1+ V), 



lh], (7) 
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where Z, 7, a = ± indicate different branches of the spec- 
trum. Here ek = (3/2)t|k|, /i = (/ii + /i2)/2 and 
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It is easy to see that, in general, h c \ < h C 2 < h c \. 

In the case of \i 7^ 0, from Eqs.([9) [TO]) one can see that 
h c i 7^ h C 2 7^ 0. In this case, one possible type of Fermi 
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FIG. 3: (color online) For symmetric doping of bilayer 
graphene, ft = 0, (a) The plot of gap solution, Aj_i, versus 
T for different levels of the mean doping, /i, at t± = 0. (b) 
The mean field phase diagram of the chiral superconductiv- 
ity showing dependence on the interlayer hopping energy, t±, 
and the temperature, T, for fi/g±i — 0,0.03,0.06,0.09 (from 
bottom to top, respectively). The line of transition from N 
to S phase is shown in which dashed and solid parts indicate 
the first and the second order transitions, respectively. 



surface topology with no effective Fermi surface for both 
I = d=l occurs when ft < h c ±. In this situation gap opens 
and pairing becomes of BCS-type (see|2k). The value of 
energy gap is 2h c i. But for \l = 0, this regime vanishes 
and the pairing becomes gapless. 

In the range h& < ft < ft C 2 different type of Fermi 
surface topology takes place. For 1=1, there is no ef- 
fective Fermi surface and excitation is gapful, whereas 
the branch I = — 1 crosses Fermi level in two points and 
consequently generates two effective Fermi surfaces, as 
depicted in Fig. [2]o. This Fermi surface topology corre- 
sponds to Sarma pairing in the phase diagram and van- 
ishes at /i = 0. 

The other type of Fermi surface topology can be real- 
ized for h C 2 < ft < h c s. In this range, there is no crossing 
point between branch I = 1 and Fermi surface, while 
the branch I = — 1 crosses Fermi surface in one point 
with a large Fermi momentum. Therefore, as is shown in 
Fig. [2k, because of large Fermi surface, the first excita- 
tion band goes to N state and the second band remains 
gapped. This topology defines new S phase which corre- 
sponds to the gapped mixed (GM) regime in the phase 
diagram. 

The last possible Fermi surface topology can be de- 
termined when h c s < ft. In this range, both branches 
Z = d=l have one effective Fermi surface (see Fig. [2]i). 
The branch 1 = —1(1) crosses Fermi surface with a large 
(small) Fermi momentum and goes to the N (S) state. 
It is interesting to note that splitting of Fermi surfaces 
by h in this case results in electron-like and holed-like 
Fermi surfaces in the absence of superconductivity. So, 
interlayer superconductivity, partially pairs n- and p-type 
momentum states. The corresponding phase will be re- 
ferred to as pn phased in the phase diagram. 

According to the above four different topologies for 
effective Fermi surfaces, the gap equation has different 



solutions. 



III. INTERLAYER SUPERCONDUCTIVITY IN 
DOPED BILAYER GRAPHENE 

We obtain gap equation by minimizing the thermody- 
namic potential 

Sls = ^±-± V ln(l + e-«<), (!2) 

with respect to the order parameters A_u, dfl/dA±i = 
0. By solving the gap equation self-consistently, one can 
obtain A_li- 

For symmetric doping of the bilayer graphene, ft = 0, 
the solution of the gap equation, A_li versus tempera- 
ture, and the phase diagram of the chiral superconduc- 
tivity in the t± — T plane, respectively, are shown in 
Fig. [3] for different levels of the mean doping. In Fig. 
[3^, the low temperature solutions of the gap equation in- 
crease with increasing the mean doping and at high mean 
doping, PA decreases monotonically with temperature re- 
sembling BCS-like behavior. This could be understood as 
follows: In the absence of doping, S state is gapless and 
the N fermions will exist as well and this weakens PA; 
in the presence of the mean doping, gap opens and PA 
strengthens. Finally, for high mean doping, Fermi sur- 
face becomes large and band structure is not important 
same as conventional s-wave symmetry, which results in 
BCS-like behavior of PA. The phase diagram of the chiral 
superconductivity in the t± — T plane is shown in Fig. 
[3b and the line of transition from N to S phases is shown 
which dashed and solid parts indicate the first- and the 
second-order phase transitions, respectively. One can see 
that low temperature region of S state increases by the 
mean doping. This is due to increase of PA by the mean 
doping [shown in Fig. [3^]. Since t± normally couples the 
two layers of the bilayer graphene, larger values of t± are 
required to vanish larger PA. 

In the case of asymmetric doped bilayer graphene, 
ft 7^ 0, the solution of the gap equation, A_li, versus 
the doping difference, ft, at zero temperature and the 
phase diagram of the chiral superconductivity in ft — T 
plane are shown in Fig. [4] for different levels of the mean 
doping and t± = 0. As depicted in Fig. 2^, at finite 
mean doping the behavior of A_li with respect to ft is 
constant up to the certain value of ft c i, decreases until 
develops a minimum at ft C 2 and then increases, at a crit- 
ical ft, PA becomes constant and at h c s abruptly goes to 
zero. These four steps of behavior of A_li correspond to 
the four different topologies of the Fermi surface as dis- 
cussed above. Approaching \i = 0, BCS and Sarma states 
disappear completely. The phase diagram of (h-T) plane 
is shown in Fig. [4}d, for different levels of the mean dop- 
ing. The stable S region is separated from N region by 
phase transition line in which dashed and solid parts of 
the lines indicate the first and the second order phase 
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FIG. 4: (color online) For asymmetric doping of bilayer 
graphene, /i 7- 0, (a) The plot of none zero solutions of 
gap equation, Aj_i, versus T for different levels of the mean 
doping fi at t± — 0. (b) The mean field phase diagram 
of the chiral superconductivity showing dependence on the 
doping difference, h, and the temperature, T, for fi/g±i = 
0, 0.03, 0.06, 0.09 (from top to bottom, respectively). The line 
of transition from N to S phase is shown in which dashed and 
solid parts indicate the first and the second order transitions, 
respectively. 
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FIG. 6: (color online) Lattice structure of trilayer graphene. 
(a) Lattice structure of ABA-stacked trilayer graphene in 
which sublattices Ai, B2 and A3 from different layers are over- 
lapped (left panel). In the right panel, interlayer correla- 
tions between electrons of sublattices A\ and B2 (B2 and A3) 
from first and second (second and third) layers are indicated 
by wavy lines, (b) Lattice structure of ABC-stacked trilayer 
graphene in which sublattices A\ and B2 from first and sec- 
ond layers and also A2 and B3 from second and third layers 
are overlapped (left panel). S correlations between interlayer 
nearest neighbor sublattices are shown in the right panel. 



IV. INTERLAYER SUPERCONDUCTIVITY IN 
TRILAYER GRAPHENE 



FIG. 5: (color online) For asymmetric doping of bilayer 
graphene, h ^ 0, phase diagram of S states in (/j, — h) plane is 
shown at zero temperature. Solid thin lines indicate bound- 
aries of different phases and thick dashed line represents first 
order phase transition from S to N state. 



transitions, respectively. At low (high) temperature and 
high (low) doping difference, the phase transition is the 
first (second) order. By increasing the mean doping, the 
range of the stable S state decreases. 

In Fig. [5j the phase diagram of S state in (/i — h) 
plane is shown at zero temperature. The solid thin lines 
indicate the boundaries of the different phases and the 
thick dashed line represents the first order phase transi- 
tion from S to N state. One can see that BCS, Sarma 
and GM phases are stable, while pn phase is unstable. 
At /i = 0, BCS and Sarma states vanish and the only 
stable state is GM phase. Upon increasing /i, the ranges 
of BCS and Sarma phases increase, with this property 
that the range of BCS state is larger than Sarma phase. 



In this section, we develop interlayer superconductivity 
to the case of trilayer graphene. Trilayer graphene con- 
sists of a monolayer graphene on the bottom of the bilayer 
graphene. There are two known manners of stacking in 
the bulk graphite which are called ABA and ABC— as 
shown in Fig. [6j In ABA stacking, the atoms of the bot- 
tommost monolayer lie exactly on the bottom of those 
of the top layer of bilayer, while for ABC stacking, one 
type of sublattice of the bottommost monolayer lies un- 
der the center of the hexagons in the bottom layer of the 
bilayer graphene. This subtle difference in stacking order 
remarkably affects on physical properties. Near the Dirac 
point, the electrons in ABC stacking behaves as massive 
fermions, while in ABA, electrons behave as both mas- 
sive and massless Dirac fermions^. We consider the total 
Hamiltonian of 7r-electrons in the trilayer graphene with- 
out S correlation^, 



H 



tri 
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ABA(C) 
_L ' 



(13) 
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FIG. 7: (color online) Solutions of gap equations Ax (1,2) ver- 
sus temperature for ABA and ABC stacking at \i — and 
t± 0. 



where H ± { ' is the nearest neighbor interlayer hop- 
ping Hamiltonian of ABA (ABC) stacked layers and has 
the following form, 

jjABA = H AB _ t± Y^(bl iita as,i,a + H.C.), (14) 

cr,i 

Hi BC = Ei B - t ± Y±(4,i,Mi,* + #- c -), (15) 

To induce the interlayer superconductivity, we introduce 
the nearest neighbor interlayer attractive interaction, 

yABA = V AB_ g±2 b\^J^ a ,a^, (16) 

a, a' ,i 

for ABA trilayer and 

yABC = yAB _ g ^ £ „ 2 , ^ , .,,/,;,;.„ , (17) 

cr,(j' ,i 

for ABC trilayer graphene. g±2 couples the electrons 
of the sublattice B 2 (A 2 ) and A 3 (B 3 ) in ABA (ABC) 
arrangement. The first term of Eqs.(pJ)J [T7|) can be de- 
coupled by Eq.(|5j). Following discussion of Sec. HU we 
determine S order parameters between the second and 
the third layers 

Af,_L2 = -0_L2(&2,a a 3,i,t + h,i,t a 3,i,l)i 

for decoupling the second term of Eq. (fl6]) and 

Ai,±2 = -g±2(a>2,i,ib 3 ^ + a 2 ,i,t&3,i,i), (19) 

for decoupling the second term of Eq. (fT7j) . In the mo- 
mentum space, the total Hamiltonian of trilayer graphene 

( H ABA(C) = H tri + yABA(C)^ ^ be diagonalized With 

the procedure of Sec. IIII1 we can calculate PA's for dif- 
ferent stacking of trilayer. Fig. [7] shows the nonzero 




FIG. 8: (color online) Finite temperature thermodynamic 
potential contours in the plane of (Axi,Ax2) showing sta- 
ble points. In the case of ABA stacking, steps (a) g±± > 
g±2 (p_Li/p_L2 = 1.02), (b) 0xi = g±2 , and (c) g ±1 < g± 2 
(g±i/g ±2=0.98) indicate first order phase transition between 
Axi and Ax2, whereas in the case of ABC stacking steps (d) 
g±i > g±2 (g±i/g ±2=3.11), (e) g ±1 = g± 2 , and (f) g ±1 < g± 2 
(g±i/g ±2=0.32) imply that second order phase transition oc- 
curs. 



stable solutions of Axi and Ax2 in terms of tempera- 
ture for pure ABA and ABC trilayer at t± = 0. Same 
as the bilayer case, there is an enhancement of PA at the 
intermediate temperature. However, the values of ABC 
trilayer order parameter are larger than those of ABA tri- 
layer, and as a result, in this case, simultaneous existence 
of Axi and Ax2 promote the interlayer superconductiv- 
ity, we will discuss below the case of ABA stacking. It 
is interesting to note that there is a mutual influence be- 
tween the interlayer PA's. The manner of effect strongly 
depends on the type of layers arrangement. This effect 
can be investigated by studying of the competition be- 
tween Axi and Ax2 with respect to the relative values 
of g±i and #X2- In Fig. [HI the contour plot of the ther- 
modynamic potential in the plane of (Axi , Ax2) for 
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ABA- and ABC-stacked layers is shown for \i = and 
t± = 0. In the ABA trilayer graphene, if g±i > g±2, the 
global (local) minimum of the thermodynamic potential 
lies on A_li (A_l2) axis (see Fig. [8fci). By increasing 
g±2 so that g±i = g±2, thermodynamic potential gives 
many stationary points where the global minimum is a 
quarter of a circle shown in Fig. [HJd. This behavior of 
global minimum is due to the symmetry of the parameter 
space. By further increasing of g±2, so that g±i < g±2, 
the minimum of the thermodynamic potential pass to 
the nonzero A^2 but zero A_li (see Fig. |5t). On the 
other hand, by varying g±i and g±2, stable point jumps 
between two axes implying first order phase transition 
between A_li and A^2- The situation for ABC-trilayer 
is different. Fig. [HJi shows that if g±i > g±2, the domi- 
nant PA is A_li and the other has a small value. On the 
other hand, the thermodynamic potential has one point 
of global minimum without local minimum. This station- 
ary point moves continually in the plane of (A_li , A^) 
by varying g±i and g±2- At the same values of g±i and 
g±2, PA's have equal magnitude (see Fig. |5fc). Finally, 
A^2 will be larger than A_li if g±i < g±2 as is shown in 
Fig. [Hp. Therefore, changing of g±s exhibits the second 
order phase transition between the order parameters. 

V. CONCLUSION 

In conclusion, we have analyzed s-wave spin-triplet su- 
perconductivity in bilayer and trilayer graphene, which 
can be realized by an interlayer pairing between chiral 
electrons of different layers. For bilayer graphene, we 
have shown that the gapless property of the interlayer S 



state of a pure sample and its associated temperature- 
induced enhancement of the S order parameter are sup- 
pressed by a finite mean doping. We have explored the 
possibility of realizing four distinct pairing regimes with 
different topologies of the Fermi surfaces by tuning the 
difference in the doping level of the two layer. These 
include the BCS regime in which two relevant branches 
of excitations are gapped and S (Fig. [2] a), the Sarma 
regime with one branch gapped and the other gapless 
both S (Fig. [2]b), gapped mixed regime with one branch 
gapped S and the other N (Fig. [2]c), and pn regime with 
both branches being gapless, one S and the other N (Fig. 
Eld). 

For trilayer graphene, we have shown that the inter- 
layer pairing strongly depends on the type of the stack- 
ing, with ABC stacking supporting higher pairing ampli- 
tude than the ABA stacking. For an imbalanced coupling 
strengths of the middle layer with the two outer layers, 
we distinct two phases which are characterized by the 
pairing between the middle layer with only one of the 
outer layers. We have found a phase transition between 
these two phases by varying the ratio of the coupling 
strengths, which is first order for the ABA stacking but 
second order for the ABC stacking. 
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